Astronomy Sz Astrophysics manuscript no. aa2318 February 2, 2008 

(DOI: will be inserted by hand later) 



On the influence of Stark broadening on Cn lines in stellar 

atmospheres 

M. S. Dimitrijevic 1 ' 2 , T. Ryabchikova 3 ' 4 , L. C. Popovic 1 ' 2 ' 5 , D. Shulyak 6 , and S. Khan 6 
1 Astronomical Observatory, Volgina 7, 11160 Belgrade 74, Serbia 

in , 

Isaac Newton Institute of Chile, Yugoslavia Branch 

3 Institute of Astronomy, Russian Academy of Science, Pyatnitskaya 48, 119017 Moscow, Russia 

4 Institute for Astronomy, University of Vienna, Tiirkenschanzstrasse 17, A- 1180 Vienna, Austria 

5 Astrophysikalisches Institut Potsdam, An der Sternwarte 16, 14482 Potsdam, Germany 
Ch ' 6 Tavrian National University, Yaltinskaya 4, 330000 Simferopol, Crimea, Ukraine 

a 1 



o 



Received / Accepted 



in 
o 



Abstract. Using the semiclassical perturbation method , electron-, proton-, and ionized helium-impact line widths 
and shifts for the nine Cri spectral lines from the 4p 7 P° — 4d 7 D multiplet, have been calculated for a perturbers 
density of 10 14 cm -3 and temperatures T = 2,500 — 50,000 K. The obtained results have been used to investigate 
' the influence of Stark broadening effect in the Cr-rich Ap star /3 CrB atmosphere on line shapes of these lines. It 

, has been found that the contribution of the proton and Hell collisions to the line width and shift is significant, 

' and it is comparable and sometimes (depending of the electron temperature) even larger than electron-impact 

contribution. Moreover, not only the Stark line width, but, depending on the electron-, proton-, and Hen density 
in stellar atmosphere also the Stark shift may contribute to the blue as well as to the red asymmetry of the same 
line. The obtained results have been used to investigate the influence of Stark broadening effect on line shapes of 
Cri lines in the atmosphere of Cr-rich Ap star f3 CrB. 
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c/) , 1. Introduction in 7 Equ (Ryabchikova et al. 2002), where electron density 

■ is high enough to favour the Stark effect. The most of Cri 

The Stark broadening is the most significant pressure and CrH lines in the optical spectral region have rather 
k> 1 broadenin S mechanism for A and B stars and this ef- gmall Stark damping constants and hence, no measurable 
; feet should be taken into account in investigation, analysis Stark wingg appeared . However, Cri lines from Ap - Ad 
g . and modelling of their atmospheres. In our previous works transitions are known to have rather large Stark damp- 
(Popovic et al. 1999, Popovic et al. 2001, Dimitrijevic ing constants acc0 rding to calculations made by Kurucz 
et al. 2003) we have shown that the Stark broadening (1993) 

may change the spectral line equivalent widths by 10-45%, We present new calculations of Cr 1 lines Stark 
hence neglecting this mechanism, a significant errors in widths and shifts based Qn the semic la SS ical perturba- 
abundance determinations may be introduced. Stark data tion approach (Sahal-Brechot 1969a,b, 1974, Dimitrijevic 
become extremely important after discovering abundance & Sahal-Brechot 1984a). Method of these calculations is 
gradients in the atmospheres of magnetic peculiar (Ap) descr ibed in Section^ magnetic synthetic spectrum calcu- 
stars (Babel 1992, Ryabchikova et al. 2002, Wade et al. lationg are descri bed in Section^ Results of Stark broad- 
2003). High resolution spectra allow us to perform strati- ening ca i cu l a tions are given in SectionO An application 
fication analysis using line profiles, and strong lines with of new data to spectrum synthes i s m magnetic Ap star 
developed wings provide us with the most accurate in- ^ CrB ig considered in Section jj 
formation about distribution of the element through the 
stellar atmosphere (see Dimitrijevic et al. 2003 for Si). 

Chromium is one of the most anomalous elements in 
Ap stars. It was shown to be concentrated in the deeper at- Calculations have been performed within the semiclas- 
mospheric layers in Ap stars (3 CrB (Wade et al. 2003) and sical perturbation formalism, developed and discussed 

in detail in Sahal-Brechot (1969ab). This formalism, as 

Send offprint requests to: M. S. Dimitrijevic well as the corresponding computer code, have been 
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optimized and updated several times (see e.g. Sahal- 
Brechot, 1974; Dimitrijevic and Sahal- Brechot, 1984a, 
Dimitrijevic, 1996). 

Within this formalism, the full width of a neutral emit- 
ter isolated spectral line broadened by electron impacts 
can be expressed in terms of cross sections for elastic and 
inelastic processes as 



2A 2 f 
W if = -^n e / vf(v)dv(£2a u ,(v} 



and the corresponding line shift as 
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Here, Xi / is the wavelength of the line originating from the 
transition with the initial atomic energy level i and the fi- 
nal level /, c is the velocity of light, n e is the electron den- 
sity, f(v) is the Maxwellian velocity distribution function 
for electrons, m is the electron mass, k is the Boltzmann 
constant, T is the temperature, and p denotes the impact 
parameter of the incoming electron. The inelastic cross 
section Oyy (v) is determined according to Chapter 3 in 
Sahal-Brechot (1969b), and elastic cross section a e i ac- 
cording to Sahal-Brechot (1969a). The cut-offs, included 
in order to maintain for the unitarity of the S-matrix, 
are described in Section 1 of Chapter 3 in Sahal-Brechot 
(1969a). 

The formulae for the ion-impact broadening param- 
eters are analogous to the formulae for electron-impact 
broadening. For the colliding ions, the validity of the im- 
pact approximation has to be checked in the far wings. 

3. Line profile calculations 

Model atmosphere calculations as well as calculations 
of the absorption coefficients were made with the lo- 
cal thermodynamical equilibrium (LTE) approximation. 
Model calculations were performed with the ATLAS9 code 
(Kurucz 1993), modified by Piskunov & Kupka (2001) to 
include individual chemical composition of the star in line 
opacity calculations (see Kupka et al. 2004). 

The next step is the calculation of the outward flux at 
corresponding wavelengths points using the given model. 
For this purpose we used SynthM code written by Khan 
(2004). This code allows to calculate synthetic spectra of 
early and intermediate type of stars taking into account 
magnetic field effects and stratification of chemical ele- 
ments. 

The computational scheme is as follows. For each line 
we find the central opacity as 



7re , _x n.„ 

aune = gfife <=T_ l_ e ht) 

mc p 



(3) 



where a u is the mass absorption coefficient at frequency 
v, e is the electron charge, g is the statistical weight, fa is 



the oscillator strength for a given transition, \ is the exci- 
tation energy, n is the number density of a corresponding 
element in a given ionization stage multiplied by partition 
function, p is the density and h is the Planck constant. 
The last factor describes stimulated emission. 

Next, we compute the total damping parameter 



T ^Jrad "i" T Stark ~t~ ^neutral- 



(4) 



Here j ra d, J stark and 7„ eutro z are the radiative broad- 
ening, Stark broadening and damping parameters due 
to neutral atom collisions respectively. The values of 
Irad, JneutraU excitation energy \ an d oscillator strength 
gf were taken from the Vienna Atomic Line Database 
(VALD) (Kupka et al. 1999). In the case of neutral 
atom broadening we assumed that perturbing particles 
are atoms of neutral hydrogen and helium only (Van der 
Waals broadening). This assumption is applicable to al- 
most all types of stars due to high hydrogen and helium 
cosmic abundances. We shall discuss the competition be- 
tween the broadenings caused by the Stark effect and neu- 
tral hydrogen collisions in the atmospheres of our template 
star (3 CrB in Section 4.2. 

In order to include Stark broadening effects we added 
the approximate formulas (see Eqs. (20) and (21) in Sec. 
4.1) in the code. The Stark width and shift are 

(e) , (p) (Hell) 

1 Stark = l S tark n e + 1 Stark n V + 1 Stark n HeII, (5) 

dstark = df tark n e + df tark n p + df^ n HeII , (6) 

where n e , n p and riHeii are the corresponding densities of 
electrons, protons and Hen ions respectively. The result- 
ing opacity profile is given by the Voigt function (Doppler 
+ pressure broadening). 

The next step is to solve transfer equation with new 
Stark damping parameters. In presence of magnetic field 
atomic level k determined by quantum numbers Jfe, L^, 
Sk splits into 2 Jk + 1 states with Mk — —Jk, ■ ■ ■ , +Jk- The 
absolute value of splitting is defined by field modulus \B\ 
and Lande factor g^ which in the case of LS coupling is 
calculated as 

l)-L fc (L fc + l) 



_ 3 Sk(Sk 

9k ~2 + 2J k (J k + l) 



(7) 



The possible transitions between splitted levels (u - upper 
and I - lower level) are allowed by selection rules 



AM = M„ 



My 




(8) 



For a normal Zeeman triplet the subscript p corresponds 
to the unshifted w component, whereas b and r correspond 
to the blue- and red-shifted a components respectively. In 
the general case of an anomalous splitting indices p, b, r 
refer to the series of the n and a components. 

The wavelength shift of the component relative to the 
laboratory line centre Aq is defined by 



AA 



eA ol-B| , 



(9) 
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The relative strengths of components S Pt r,b in accordance 
with Sobelman (1977) are proportional to 



1 Jl 



M„. M„ - Mi Mi 



(10) 



where the last structure is a 3j-symbol. The normalization 
of components is done so that 



(11) 



In general case the polarized radiation can be described 
by means of Stokes IQUV parameters. The transfer equa- 
tion is 



dl 

dz 



KI+J. 



Here J = (I, Q, U, V) T is the Stokes vector, K is the ab- 
sorption matrix and J is the emission vector, 

K = CC C 1 + ^ ai ine & n ne , (13) 
lines 

J v = a c B v (T)e + B V (T) ^ a Une $> Une e , (14) 

lines 

where 1 is the 4x4 unit matrix, eo = (1,0, 0,0) T , a c 
is the continuum absorption coefficient, au ne is the line 
central opacity for zero damping and zero magnetic field 
given by equation (3). Here we also assume LTE so the 
source function is equal to Planck function B U (T). 

The line absorption matrix *&u ne is created using ab- 
sorption profiles 4>j and anomalous dispersion profiles tpj 
of 7r and a± components in accordance with Rees et al. 
(1989). These profiles are calculated as (for j — p, b, r) 



N, 



E 5 ^ 



AAj 

AAr 



A', 



in — 1 



(15) 



(16) 



where H(a, u) and F(a, u) are the Voigt and Faraday- 
Voigt functions, 



7A§ 



47TcAAl3 ' 

A — Ao + dstark 

AAd ' 



(17) 



(18) 



The Stark shift dg tar ^ and the damping parameter 7 have 
been found from (6) and (4), respectively. The Doppler 
width AAo is 



A 2kT 



A\ n = -J—+g 7 
c V rriA 



(19) 



where vtla is the mass of the absorber atom and £t is the 
microturbulent velocity. 



The transfer equation 1)1 2|l is solved by DELO method 
(Rees et al. 1989) with quadratic approximation for the 
source function which was suggested by Socas-Navarrro 
et al. (2000) and written in Fortran by Piskunov & 
Kochukhov (2000). Also, we used fast Humlicek (1982) 
algorithm for approximations of Voigt and Faraday- Voigt 
functions. 

The convolution of the synthetic spectra due to rota- 
tion of the star was produced by IDL procedure of Valenti 
& Anderson implemented in the rotate IDL code by 
N. Piskunov. 

4. Results 



(12) 4.1. Stark broadening data 



The atomic energy levels needed for Stark broadening cal- 
culations were taken from Wiese and Musgrove (1989). 
Oscillator strengths have been calculated by using the 
method of Bates and Damgaard (1949) and tables of 
Oertel and Shomo (1968). For higher levels, the method 
described in van Regemorter, Binh Dy and Prud'homme 
(1979) has been applied. 

The spectrum of neutral chromium is complex and not 
known well enough for a good calculation of the considered 
lines. The principal problem is the absence of the reliable 
experimental data for / levels. However, by inspecting the- 
oretically predicted Cr I levels in Moore (1971), one can 
see that the 4/ level is missing. In accordance with the 
decrease of distance between existing s,p, and d levels, 
we estimated that the eventual 5/ level should be around 
52000 cm -1 and not closer than 50000 cm -1 . We checked 
results without 5/ level and with a Active 5/ level at 52000 
and 50000 cm -1 . In all cases, line widths differed by less 
than 1%, while the shift varies within the limits of several 
percents. Consequently we performed calculations without 
the contribution of / energy levels. Since the average esti- 
mated error of the semiclassical method is ±30%, due to 
additional approximations and uncertainties, we estimate 
the error bars of our results to be ±50%. 

Our results for electron-, proton-, and ionized helium- 
impact line widths and shifts for the nine considered Cri 
spectral lines for a perturbed density of 10 14 cm~ 3 and 
temperatures T = 2,500 — 50,000 K, are shown in Table 
1. For perturber densities lower than those tabulated here, 
Stark broadening parameters vary linearly with perturber 
density. The nonlinear behavior of Stark broadening pa- 
rameters at higher densities is the consequence of the in- 
fluence of Debye shielding and has been analyzed in detail 
in Dimitrijevic and Sahal-Brechot (1984b). 
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Table 1. Stark broadening parameters for Cn 4p - 4d spectral lines. This table shows electron-, proton-, and ionized 
helium-impact broadening parameters for Cn for a perturber density of 10 14 cm" 3 and temperatures from 2,500 up to 
50,000 K. The quantity C (given in A cm~ 3 ), when divided by the corresponding full width at half maximum, gives an 
estimate for the maximum perturber density for which tabulated data may be used. For higher densities, the isolated 
line approximation used in the calculations breaks down. The asterisk identifies cases for which the collision volume 
multiplied by the perturber density (the condition for validity of the impact approximation) lies between 0.1 and 0.5. 
When this value is larger than 0.5, the corresponding results have been omitted. FWHM(A) denotes full line width at 
half maximum in A, while SHIFT(A) denotes line shift in A. 



PERTURBERS ARE: 


ELECTRONS 


PROTONS 


HELIUM IONS 


TRANSITION 


T(K) 


FWHM(A) 


SHIFT(A) 


FWHM(A) 


SHIFT(A) 


FWHM(A) 


SHIFT(A) 


Cr I 


■'Pg - 'D x 


2500. 


0.890E-02 


-0.205E-02 


0.461E-02 


-0.379E-02 


0.344E-02 


-0.299E-02 






5000. 


0.772E-02 


-0.146E-02 


0.551E-02 


-0.437E-02 


0.388E-02 


-0.346E-02 




5276.07 A 


10000. 


0.655E-02 


-0.104E-02 


0.678E-02 


-0.505E-02 


0.439E-02 


-0.396E-02 


C = 


0.43E+15 


20000. 


0.548E-02 


-0.767E-03 


0.832E-02 


-0.592E-02 


0.503E-02 


-0.451E-02 






30000. 


0.493E-02 


-0.660E-03 


0.915E-02 


-0.651E-02 


0.547E-02 


-0.488E-02 






50000. 


0.431E-02 


-0.559E-03 


0.993E-02 


-0.719E-02 


0.602E-02 


-0.544E-02 


Cr I 


'Pa - 'D 2 


2500. 


0.193E-01 


-0.216E-02 














5000. 


0.164E-01 


-0.155E-02 












5275.75 A 


10000. 


0.136E-01 


-0.108E-02 










C = 


0.11E+14 


20000. 


0.111E-01 


-0.790E-03 


*0.717E-01 


-0.298E-01 










30000. 


0.983E-02 


-0.675E-03 


*0.702E-01 


-0.273E-01 










50000. 


0.841E-02 


-0.569E-03 


*0.664E-01 


-0.239E-01 


*0.135E-01 


-0.346E-01 


Cr I 


''PS - 'D 3 


2500. 


0.240E-01 


-0.939E-03 


*0.216E-01 


-0.138E-01 


*0.131E-01 


-0.102E-01 






5000. 


0.210E-01 


-0.677E-03 


*0.253E-01 


-0.170E-01 


*0.151E-01 


-0.124E-01 




5275.28 A 


10000. 


0.178E-01 


-0.473E-03 


*0.274E-01 


-0.193E-01 


*0.169E-01 


-0.150E-01 


C = 


0.60E+14 


20000. 


0.147E-01 


-0.401E-03 


0.280E-01 


-0.193E-01 


*0.170E-01 


-0.179E-01 






30000. 


0.131E-01 


-0.393E-03 


0.283E-01 


-0.180E-01 


*0.157E-01 


-0.193E-01 






50000. 


0.112E-01 


-0.393E-03 


0.289E-01 


-0.155E-01 


*0.127E-01 


-0.200E-01 


Cr I 


Y P§ - 'D a 


2500. 


0.195E-01 


-0.219E-02 














5000. 


0.165E-01 


-0.156E-02 












5298.49 A 


10000. 


0.137E-01 


-0.109E-02 










C = 


0.12E+14 


20000. 


0.112E-01 


-0.801E-03 


*0.723E-01 


-0.301E-01 










30000. 


0.992E-02 


-0.685E-03 


*0.708E-01 


-0.276E-01 










50000. 


0.849E-02 


-0.578E-03 


*0.670E-01 


-0.242E-01 


*0.136E-01 


-0.349E-01 


Cr I 


'P°i - 'D 3 


2500. 


0.243E-01 


-0.950E-03 


*0.218E-01 


-0.140E-01 


*0.133E-01 


-0.103E-01 






5000. 


0.212E-01 


-0.687E-03 


*0.256E-01 


-0.171E-01 


*0.153E-01 


-0.125E-01 




5298.02 A 


10000. 


0.179E-01 


-0.481E-03 


*0.276E-01 


-0.195E-01 


*0.171E-01 


-0.151E-01 


c = 


0.60E+14 


20000. 


0.148E-01 


-0.408E-03 


0.283E-01 


-0.195E-01 


*0.171E-01 


-0.180E-01 






30000. 


0.132E-01 


-0.400E-03 


0.286E-01 


-0.182E-01 


*0.158E-01 


-0.195E-01 






50000. 


0.113E-01 


-0.400E-03 


0.291E-01 


-0.157E-01 


*0.128E-01 


-0.201E-01 



After testing the density dependence of Stark param- — [A] = C2 • (A + C ■ T B 
eters, we have found that the width and shift are linear n e 
functions of density for perturber densities smaller than 
10 16 cm~ 3 and can be scaled by the simple formula: 



[W 1 d) N = (W 1 d\( — ) 1 (20) 

where (W, d)jv are the width and shift at a perturber den- 
sity N (cm -3 ), and (W, d)o are width and shift given in 
Table 1, respectively. 

In order to simplify the use of Stark broadening data 
in the codes for stellar spectral synthesis, we have found 
an analytical expression for Stark widths and shifts 

-[A] = Cl .(A + T s ), (21) 



The constants c\, c 2 , A, B and C are given in Table 2. 
We take T as T/(10000if). The analytical fit is present 
in Figs. 1 and 2. As one can see from the Figs. 1 and 2, 
the analytical assumption good fit the calculated value in 
both, Stark widths and shifts. Here we should note that in 
some lines (as e.g. in Cr I A = 5329.14 A) the analytical fit 
in the case of proton and He n impact broadening param- 
eters not exactly fit the calculated values (see e.g. Fig. 2, 
fit for proton- impact shift), but the differences are smaller 
than the expected error of our calculations (« ±50%). 
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PERTURBERS ARE: 


ELECTRONS 


PROTONS 


HELIUM IONS 


TRANSITION 




FWHM(A) 


SHIFT(A) 


FWHM(A) 


SHIFT(A) 


FWHM(A) 


SHIFT(A) 


Cr I 


Y PS - Y D„ 

r 3 ^4 


2500. 


0.182E-01 


0.452E-03 


0.109E-01 


-0.737E-02 


*0.683E-02 


-0.564E-02 






5000. 


0.156E-01 


0.277E-03 


0.130E-01 


-0.902E-02 


*0.788E-02 


-0.664E-02 






10000. 


0.131E-01 


0.129E-03 


0.144E-01 


-0.102E-01 


0.909E-02 


-0.787E-02 


c = 


0.12E+15 


20000 


108E-01 


-0.434E-04 


168E-01 


-0.945E-02 


956E-02 


-0.934E-02 






30000. 


0.956E-02 


-0.141E-03 


0.172E-01 


-0.830E-02 


0.945E-02 


-0.100E-01 






50000. 


0.822E-02 


-0.240E-03 


0.194E-01 


-0.603E-02 


0.793E-02 


-0.100E-01 


Cr T 




2500. 


0.245E-01 


-0.956E-03 


*0.220E-01 


-0.141E-01 


*0.134E-01 


-0.104E-01 






5000. 


0.214E-01 


-0 701E-03 


*0.259E-01 


-0 173E-01 


*0 155E-01 


-0.127E-01 






10000. 


0.181E-01 


-0.492E-03 


*0.280E-01 


-0.197E-01 


*0.173E-01 


-0.153E-01 


c = 


0.61E+14 


20000. 


0.150E-01 


-0.419E-03 


0.286E-01 


-0.197E-01 


*0.173E-01 


-0.182E-01 






30000. 


0.134E-01 


-0.411E-03 


0.289E-01 


-0.184E-01 


*0.160E-01 


-0.197E-01 






50000. 


0.115E-01 


-0.410E-03 


0.295E-01 


-0.159E-01 


*0.130E-01 


-0.204E-01 


Cr T 




2500. 


0.184E-01 


460E-03 


0.110E-01 


-0 746E-02 


*0.692E-02 


-0 570E-02 






5000. 


158E-01 


276E-03 


0.132E-01 


-0.911E-02 


*0 798E-02 


-0.672E-02 




ci^oq 14 A 

<J i_/ • _1_ t: A 


10000. 


0.132E-01 


0.127E-03 


0.146E-01 


-0.103E-01 


0.920E-02 


-0.796E-02 


c = 


0.12E+15 


20000. 


0.109E-01 


-0.484E-04 


0.170E-01 


-0.956E-02 


0.967E-02 


-0.946E-02 






30000. 


0.968E-02 


-0.147E-03 


0.174E-01 


-0.840E-02 


0.956E-02 


-0.101E-01 






50000. 


0.832E-02 


-0.248E-03 


0.197E-01 


-0.610E-02 


0.802E-02 


-0.101E-01 


Cr I 


'PI - 'D 5 


2500. 


0.920E-02 


0.162E-02 


0.499E-02 


0.407E-02 


0.370E-02 


0.321E-02 






5000. 


0.800E-02 


0.979E-03 


0.599E-02 


0.470E-02 


0.417E-02 


0.372E-02 




5328.32 A 


10000. 


0.678E-02 


0.474E-03 


0.739E-02 


0.544E-02 


0.473E-02 


0.426E-02 


C = 


0.40E+15 


20000. 


0.567E-02 


0.103E-03 


0.902E-02 


0.641E-02 


0.543E-02 


0.485E-02 






30000. 


0.509E-02 


-0.604E-04 


0.987E-02 


0.709E-02 


0.590E-02 


0.526E-02 






50000. 


0.445E-02 


-0.208E-03 


0.106E-01 


0.784E-02 


0.647E-02 


0.589E-02 



0.25 




0.05 - 



1 


2 3 
T(in 10000 K) 


4 


5 








— — -4 


- ,4' 








-*' . ~e — 


• 


2 3 


4 


5 



T (in 10000 K) 

Fig. 1. The analytic fit of Cr I A = 5329.14 A (upper Fig.) 
and A = 5328.32 A (down) Stark widths due to impact 
with: electrons (open circles), protons (full circles) and 
Hen (full triangles). 

4.2. The Stark broadening effect on the shape of Cn 
lines 

In spite of the rather large Stark damping constants the ef- 
fect is not observable in stars with solar Cr abundance. In 
hot stars where electron and proton densities are high, Cn 
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Fig. 2. The same as in Fig. 1, but for Stark shift. 

lines considered are generally very weak, while in cooler 
stars (solar type) where these lines are strong enough other 
broadening effects are more significant. The only chance 
to look at Stark effect is the stratified atmosphere of Cr- 
rich Ap star. The best example is the well known magnetic 
star (3 CrB. 

Our spectral analysis was based on the spectrum of 
CrB obtained in February 1998 with MuSiCoS spec- 
tropolarimetcr mounted on the 2 m telescope at Pic du 
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Table 2. The parameters A, B and C of the approximate formulae for Stark widths and shifts. 



Line 


5276.07 A 


5275.75 A 


5275.28 A 


5298.49 A 


5298.02 A 


5297.38 A 


5329.78 A 


5329.14 A 


5328.32 A 


WIDTH 


electrons 


Cl 


IE- 16 


IE- 15 


IE- 15 


IE- 15 


IE- 15 


IE- 15 


IE- 15 


IE- 15 


IE- 16 


A 


-0.3474 


-0.86156 


-0.8218 


-0.86039 


-0.82019 


-0.8674 


-0.81818 


-0.86596 


-0.325497 


B 


-0.15596 


-0.03684 


-0.04346 


-0.03716 


-0.04412 


-0.0337 


-0.04412 


-0.03410 


-0.161576 


protons 


Cl 


IE- 16 


IE- 15 


IE- 15 


IE- 15 


IE- 15 


1E-15 


IE- 15 


IE- 15 


IE- 15 


A 


-0.3194 


-0.23723* 


-0.7393 


-0.2312* 


-0.7368 


-0.85357 


-0.7338 


-0.85169 


-0.924767 


B 


0.1789 


-0.06276* 


0.02231 


-0.062768* 


0.02238 


0.02723 


0.02281 


0.027692 


0.01959 


Hen 


Cl 


IE- 16 




IE- 15 




IE- 15 


IE- 16 


IE- 15 


IE- 16 


IE- 16 


A 


-0.5525 




-0.84377 




-0.84219 


-0.14236 


-0.84037* 


-0.13199 


-0.51847 


B 


0.08565 




0.01192 




0.011451 


0.11010 


0.01177* 


0.110894 


0.092283 


SHIFT 


electrons 


C2 


IE- 16 


IE- 16 


IE- 17 


IE- 16 


IE- 17 


IE- 17 


IE- 17 


IE- 17 


IE- 16 


A 


-3.4929 


-6.7617 


-81.5107 


-2.5259 


-9.0041 


49.6736 


-16.5511 


40.5533 


23.3420 


B 


0.01451 


0.00793 


0.00221 


0.02206 


0.02144 


0.00526 


0.01131 


0.00586 


0.00261 


C 


3.3728 


6.6362 


80.9250 


2.3985 


8.4079 


-49.60704 


15.9454 


-40.4305 


-23.2800 


protons 


C2 


IE- 16 


IE- 15 


1E-15 


IE- 15 


IE- 15 


IE- 16 


IE- 15 


1E-16 


1E-16 


A 


-0.0170 


0.25969* 


3.0908 


-3.3759* 


0.5500 


-2.4538 


3.7646 


18.8373 


50.1882 


B 


0.22391 


-0.122027* 


0.00215 


0.00846* 


0.00990 


0.02228 


0.00188 


-0.00173 


-0.002557 


C 


-0.4916 


-0.6076* 


-3.2607 


3.0621* 


-0.7217 


1.6060 


-3.9381 


-19.6934 


-49.62339 


Hen 


C2 


IE- 16 




1E-15 




IE- 15 


1E-16 


IE- 15 


1E-16 


1E-16 


A 


-0.0174 




1.7973 




2.4949 


37.4669 


2.0440 


36.5920 


0.0405 


B 


0.20567 




0.01778 




0.01314 


0.00418 


0.01604 


0.00431 


0.22021 


C 


-0.3771 




-1.9474 




-2.6462 


-38.2561 


-2.1971 


-37.3904 


0.3832 



Midi observatory (R=35000). This spectrum was used by 
Wade et al. (2001,2003) in stratification analysis. Starting 
Cr and Fe distributions were taken from this analysis 
and were slightly corrected using a set of spectral lines 
of Cri Cr II Fei Fe II with different excitation energies of 
the lower level. The adopted distributions are shown in 
Fig. |31 Fe stratification is necessary to account correctly 
for blends. 

Synthetic spectrum calculations in the regions of Cri 
lines considered in the present paper was performed fol- 
lowing the method described in Section |3| The surface 
magnetic field B s =5.4 kG was derived from the magneti- 
cally splitted lines and used in all calculations. In the pro- 
cess of calculations we found that Cri lines of our interest 
have incorrect wavelengths with the shifts up to 0.05 A. 
We adjusted wavelengths calculating solar synthetic spec- 
trum and comparing it with the Solar Flux Atlas (Kurucz 
et al. 1984). Improved wavelengths of Cri lines derived by 
this way agree within 0.005 A and better with the precise 
measurements by J. E. Murray (1992), kindly provided 
to us by R. Kurucz 1 . We placed Murray's wavelengths in 
Table 1 and 2 and used them in spectral synthesis. We also 
estimated the contribution of Van der Waals broadening 

1 http:/ /cfaku5. cfa.harvard.edu/atoms/2400/ 




log ^ 



Fig. 3. Empirically derived Cr and Fe distributions in the 
atmosphere of magnetic Ap star (5 CrB. 



and found it to be by the order of magnitude smaller then 
Stark broadening through the whole atmosphere above 
^°ff( T 500o) = —0.2 with the rapid decrease below the pho- 
tosphere. 

We varied Stark broadening parameters given in 
Tabled to achieve the best fit of the calculated Cri line 
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Fig. 4. A comparison between synthetic spectrum calcula- 
tions and the observed spectrum (thick line) of magnetic 
Ap star (3 CrB in the regions of Cri lines 5297 A (top 
panel). 5276 A (middle panel) and 5329 A (bottom panel). 
Full thin line - calculations with stratification shown in 
Fig |2land Stark broadening data from Table H decreased 
by 70 %; dashed line - the same Stark broadening data but 
homogeneous Cr and Fe abundances (see text); dotted line 
- stratified Cr and Fe abundances, but Stark broadening 
calculated with the approximation formula. 



profiles to the observations. Fig. 0] shows a comparison 
between observations and calculations in the region of 
3 groups of Cri lines. For comparison we present syn- 



thetic calculations with homogeneous mean abundances 
log(CY/iV tot ) = -4.80 and log(Fe/N to t) = -4.15 which 
were used in model atmosphere calculations (Kupka et 
al. 2004). To fit Cri line wings we need to decrease cur- 
rent Stark widths by 60-70 %. The same order of overes- 
timate of Stark damping constants, calculated within the 
semiclassical perturbation formalism, was obtained for Sii 
lines (see Dimitrijevic et al. 2003). Any change in Cr dis- 
tribution, which may be responsible for line profiles, too, 
is critical for nearby Cri A 5296.691 A, which is insensi- 
tive to Stark broadening mechanism due to small value of 
Stark damping constant and is sensitive only to Cr dis- 
tribution in stellar atmosphere. We also made a compar- 
ison with synthetic spectrum calculated using stratified 
Cr distribution and an approximation formula for Stark 
broadening (Cowley 1971). This formula is used in the 
cases when Stark broadening data, experimental or the- 
oretical, are not available. Evidently, Stark broadening is 
overestimated when calculations are performed with this 
approximation formula. 

For detailed investigation of the influence of Stark 
shifts on the observed line profile we calculated theoretical 
spectra with and without shifts for two separate Cri lines 
5297.377 A and 5298.494 A in homogeneous and stratified 
atmosphere of the non-rotating star. These two lines have 
the smallest and largest shift values, therefore one can 
make a relative comparison. Fig. [5] illustrates an influence 
of the Stark shift on line shape. The total shift in line po- 
sition due to the Stark broadening effect is negligible for 
Cri 5297.377 A line both in homogeneous and stratified 
atmosphere. For the second line Stark broadening leads 
to theoretically measurable shift in the line position, in 
particular, in the case of Cr stratification resulting in a 
small line asymmetry. Unfortunately, in spectra of real 
stars effect is difficult to measure due to other broadening 
mechanisms (rotation, magnetic field) and mainly due to 
severe blending in the regions with Cri lines from Ap — Ad 
transitions. 

5. Conclusions 

Stark broadening parameters for Cr i spectral lines from 
4p7p0 _ 4^7 _£) multiplet have been calculated and the in- 
fluence of Stark broadening effect in stellar atmosphere for 
these lines, has been investigated. From our investigation 
we can conclude: 

(i) Although they belong to a same multiplet, the 
widths and shifts of the different lines can be quite dif- 
ferent. 

(ii) The contribution of the proton and He II collisions 
to the line width and shift is significant, and it is compa- 
rable and sometimes (depending of the electron tempera- 
ture) even larger than electron-impact contribution. 

(iii) Depending on the electron-, proton-, and Hen 
density in stellar atmosphere the Stark shift may con- 
tribute to the blue as well as to the red asymmetry of 
the same line (see Fig. 2). 



A/r q n;^,;+^;; Q ,r;^ o+ „i . Stark broadening of Cn lines 



1 1 1 

CrI 


1 i 

5297377 Cr I 


1 

5298.494 




f \ 


r 


i 


J 

1,1, 


i 



!96.5 5297 5297.5 5298 5298.5 5299 5299.5 



1 i 1 i 1 i 1 i 1 r 

Cr 1 5297.377 Cr I 5298.494 




Wavelength, A 



Fig. 5. Influence of Stark shift on Cn line shapes in homo- 
geneous (upper panel) and in stratified (lower panel) at- 
mosphere. Calculations with both Stark widths and shifts 
are shown by solid lines while those without Stark shifts 
are shown by dashed line. 

(iv) To fit well Cr 1 line wings we need to decrease 
the calculated Stark widths by 60-70%, which is the same 
order of overestimation as for Si 1 lines (Dimitrijevic et al. 
2003). The approximation formula of Cowley (1971), used 
in the cases where the adequate semiclassical calculation 
is not possible due to the lack of reliable atomic data, 
predicts also overestimated influence of Stark broadening 
in comparison with observations. 
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